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Abstract 

The first three coefficients of the Taylor's series expansion of the vevtor pion form 
factor as a function of the momentum transfer are evaluated using the experimental 
data on the pion form factor and the P-wave tttt phase shifts. The real part of 
the form factor as a function of energy is also calculated by dispersion relation. 
Comparisons there results with Chiral Perturbation Theory and unitarized models 
are given. 
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The purpose of this note is three-fold. First, a systematic procedure is given to calcu- 
late the coefficients of the Taylor's series expansion of the form factor around the origin 
s = in terms of experimental data. Second the real parts of the form factor for a wider 
range of energy are calculated using also experimental data. Third, these results are 
compared with those given by existing theoretical calculations in order to evaluate their 
reliability. It turns out that only models satisfying the elastic unitarity, giving rise to the 
p resonance, are valid. Chiral Perturbation calculations at one and two-loop disagree, at 
very low energy, with the Taylor's series expansion and, at moderate energy, with the 
dispersion relation results. 

Chiral Perturbation Theory (ChPT) |], @, 0, f| is a well-defined perturbative procedure 
allowing one to calculate systematically low energy phenomenon involving soft pions. It 
is now widely used to analyze the low energy pion physics not only when the interaction 
is not strong but also in the presence of the resonance as long as the energy region of 
interest is sufficiently far from the resonance. In this scheme, the unitarity relation is 
satisfied perturbatively order by order. 

In this note we want to examine critically the ChPT approach and wish to emphasize 
that unless that this calculation was followed by a unitarisation procedure, its results for 
a number of physical processes would not be meaningful. This is so because in any pertur- 
bative scheme, it is important that the magnitude of the calculated term should be much 
larger than that of the same order which is uncalculable because of its non-perturbative 
origin. To make this point clear, let us examine the one-loop ChPT calculation of the 
vector pion form factor V(s) |]: 

V perL (s) = 1 + - + wAr 2 ((s - ±rnl)H™{s) + ^) (1) 
s R 96% 2 f£ 3 

where f n = 0.93GeV and the r.m.s. radius of the vector form factor is related to sr by 
the definition V (0) = |ry = 1/sr. The function H nn (s) is given by: 



/ o _ 4 m 2 y/s + Js — Ami I s — Am 2 

fors > Am 2 ; for other values of s, H niT {s) can be obtained by analytic continuation. 

From this definition, it is clear that the third term on the r.h.s. of Eq. (JTJ), being an 
analytic function with the same singularity as V(s), can be expanded as a Taylor's series 
at s = with a leading term behaving as s 2 at small s. For small s it is a reasonable 
approximation to take only the first few terms of this series. 

The unitarised version of Eq. ([!]), obtained by the inverse amplitude, the Pade ap- 
proximant or the N/D methods, is given by [[7]. |J: 

V(S) = 1 - s/sr - ^{(s - Aml)H W7r (s) + 2s/3} (3) 

It is obvious that Eq. (||) has the Breit-Wigner resonance character while that from 
Eq. ([!]) does not, although their amplitude and first derivative are identical at s — 0. 
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Furthermore, if the parameter sr was fixed by the the r.m.s. radius, the p mass would 
come out to be slightly low compared measured p mass. Neglecting this last problem at 
the moment, the Taylor's series expansion around s = reveals that Eq. (0) gives rise to 
a coefficient of the s 2 term as (1/sr) 2 ~ A.OGeV~ A which is much larger than that coming 
from the third term of Eq. ([I]), l/(9607r 2 m 2 / 2 ) ~ O.Q3GeV~ 4 . This is the signal of the 
failure of the perturbation method. While it is difficult to detect the presence of this term 
using experimental data at low energy, this failure would show up at higher energy and 
with more accurate data. 

It could be argued, however, that this discrepancy could come from the unitarisation 
scheme, because Eq. (|3]) fits the data only approximatively. It is therefore desirable to 
have a direct proof that this discrepancy is real and only depends on the experimental 
data. 

Because the vector pion form factor V(s) is an analytic function with a cut from 4m 2 
to oo, the n th subtracted dispersion relation for V(s) reads: 



„_i i s n f°° ImV(z)dz 



V(s)=a + a 1 s + ...a n _ 1 s n L + — — v , (4) 

7T Jiml Z n (z - S—ie) 

where n > and, for simplicity, the series around the origin is considered. The polynomial 
on the R.H.S. of this equation will be referred in the following as the subtraction constants 
and the last term on the R.H.S. as the dispersion integral (DI). The evaluation of DI as 
a funtion of s will be done later. Notice that a n = V n (0)/n\ is the coefficient of the 
Taylor series expansion for V(s), where V n (0) is the nth derivative of V(s) evaluated at 
the origin. The condition for Eq. (|]) to be valid was that, on the real positive s axis, the 
limit s~ n V(s) — > as s — > oo. By the Phragmen Lindeloff theorem, this limit would also 
be true in any direction in the complex s-plane and hence it is straightforward to prove 
Eq. (H). The coefficient a n+m of the Taylor's series is given by: 

1 r°° ImV(z)dz 

7T J Ami ^ n+m+1 ) ^ ' 



"n+m 



where m > 0. The meaning of this equation is clear: under the above stated assumption, 
not only the coefficient a n can be calculated but all other coefficients a n+m can also be 
calculated. 

Below the inelastic threshold, from unitarity of the S-matrix, V(s) must have the 
phase of the P-wave elastic tttt phase shifts From the available experimental data, the 
inelastic effect will manifest only at an energy above 1.3GeV [fTPf . Therefore, below this 
energy one has: 

ImV(z) =\V(z) \ sin 5 (z) (6) 

where 6 is the strong elastic P-wave phase shifts. (There is an ambiguity of ± sign 
on the RHS of Eq. the choice of the sign + was made is due to the normalisation 
V(0) = 1 and the phenomenology, e.g. the rms pion radius is positive). Because the real 
and imaginary parts are related by dispersion relation, it is important to know accurately 
ImV(z) over a large energy region. Below 1.3 GeV, ImV(z) can be determined accurately 
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because the modulus of the vector form factor [p| , |T2|j and the corresponding P-wave irir 
phase shifts are well measured [13|, [14], |15 |. 



Using the experimental data on the pion form factor and the corresponding nn phase 
shifts, ImV(z) with experimental errors is given in Fig. 1. Similarly in Fig. 2, ReV(s) 
with experimental errors is given using the expression ReV(s) =\ V(s) | cos S. 

One first shows how the coefficients of the Taylor's series around s = can be evaluated 
in terms of ImV(s) then later, the experimental ImV(s) and ReV(s) will be used to test 
the validity of various theoretical model calculations . 

Following the usual definition, V(s) = l + ^<r v >s + cs 2 + ds 3 + ... one has 
< ry >= 6ai,c = ci2,d = as etc. If one makes a weak assumption that V(s)/s — > as 
s — > oo, using Eq. (|5|) one gets the following results: 

< r 2 v >= 0.44 ± 0.015/m 2 ; c = 3.90 ± OAOGeV' 4 ; d = 9.70 ± 0.40GW" 6 (7) 

where the upper limit of the integration is taken to be 1.7GeV 2 . From the 2 ir threshold 
to 0.5GeV the experimental data on the the phase shifts are either poor or unavailable, 
an extrapolation procedure has to be used. For this purpose the results of models 1 and 
2, to be discussed later, were used for ImV(s). They contribute, rspectively, 5%, 15% 
and 30% to the a%, ai and a 3 sum rules. 

Because of the assumed high energy behavior, one cannot calculate here the pion 
charge ao- From these results the radius of convergence for the Taylor's series is (much) 
less than lGeV 2 . In fact, a reasonable approximation for this series is the Taylor's series 
expansion of V(s) ~ 1/(1 — ais). This result is not surprising because it is the zero width 
approximation for Eq. (|3|). Away from the p resonance, it is a better approximation than 
the ChPT calculations for the vector pion form factor. 

If one was willing to make a stronger assumption that the form factor vanished asymp- 
totically, then one would be able to calculate the pion charge ao and of course also higher 
derivatives. A straightforward calculation gives ao = 1.02 ± .08, where the upper limit of 
the integration is taken to be 1.7GeV 2 . The Ward identity requires a to be exactly unity. 
This calculation is to illustrate the fact that the p resonance gives the major contribution 
and almost saturates the sum rule for the charge; it will not be used in the following. The 
higher derivative sum rules are much less sensitive to the high energy behavior of ImV(z) 
because of the corresponding weight factor. 

It is possible to estimate the high energy contribution to the sum rules by fitting the 
asymptotic value of the pion form factor by the expression, V(s) = — ln(— s/s p ) which 
fits equally well the large energy behavior of the time-like and space-like pion form factor, 
then the integration of the sum rule from 1.7GeV 2 to infinity can be estimated. The 
r.m.s. radius is then increased by 2% and is completely negligible for the values of c and 
d. (The additional contribution to the a sum rule is 15%). 

Notice that the determination of the pion r.m.s. agrees well with its direct experimen- 
tal value of 0.439 ± .008/m 2 |L6] and with comparable errors. The determination of the 
values of c given by emperical fits to the data using the vector meson dominance models, 
quoted as having a large error of 30%, is due to the use of models having different values 
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for the form factor at the p mass || 12 1. The present determination of c as well as of d 
has an error of only a few percents. 

It is now possible to compare these experimental results with existing calculations. 

The first model is obtained using the inverse amplitude method for the pion form 
factor at one-loop level. It can also be obtained by using the (0,1) Pade approximant 
method for the one-loop ChPT calculation. The result is well known: by fitting the time 
like the pion form factor using the experimental p mass, the r.m.s. radius is 10% too 
small. The calculated p width is 0.156GeV which is satisfactory but the maximum value 
of square of the modulus of the pion form factor is only 31 which is too low compared 
with the experimental data value of 40.5 ± 0.6 as given directly by the r decay data. 

This result is based on the assumption of the elastic unitarity which should be correct 
for an energy below 1.3 GeV but is certainly incorrect above this energy. One should 
phenomenologically correct for this strong assumption. To do this one can simulate the 
inelastic effect (and possibly the two-loop singularity) by the presence of the polynomial 
ambiguity in the phase representation of the form factor i.e. a zero in the form factor. 
Instead of Eq. (|3J), the time-like pion form factor data can now be fitted by |HJ : 



y, x 1 + as/s p 

{S) 1 - s/s R - ^{(s - Aml)H^{s) + 2s/3} 1 ' 

where /„. = 0.093GeV, and sr is related to the p mass squared s p by requiring that the 
real part of the denominator of Eq. (§) vanishes at the p mass. 

The experimental data can be fitted with a p mass equal to Q.17?>GeV and a = 0.14. 
This is a two parameter fit to the experimental data and there are three predictions: the 
pion r.m.s. radius, the p width and the value of the form factor at the p mass (or the p 
leptonic width). These results are in excellent agreement with the data ]T2, |Tj|: the p 



width, defined as the derivative of the phase shift at the p mass, is equal to 0.156GeV, 
the rms radius is predicted to be 0.45 ± 0.01/m 2 and the maximum value of the modulus 
squared of the pion form factor at the p resonance is 39.2. 

The second model, which is more complicated, but is more complete because it is 
based on the two- loop ChPT calculation with unitarity taken into account. It has the 
singularity associated with the two loop graphs. By using the same inverse amplitude 
method as was done with the one-loop amplitude, but by generalizing this method to 
two-loop calculation, Hannah has recently obtained a remarkable fit to the pion form 
factor in the time-like and space-like regions. His result is equivalent to the (0,2) Pade 



approximant method as applied to the two- loop ChPT calculation |L8]|. Unlike Eq. flg) 
which is a two parameter fit, Hannah calculation is a three parameter calculation which 
gives the correct maximum value of the pion form factor, the p mass and width and the 
pion r.m.s. radius. 

As can be seen from Figs. 1 and 2, the imaginary and real parts of these two models 
are very much in agreement with the data. A small deviation of ImV(s) above 0.9GeV 
is due to a small deviation of the phases of V(s) in these two models from the data of the 
P-wave 7T7T phase shifts. They both give also correct results for c and d as given by the 
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sum rules, Eq.(^). For < s < OAGeV 2 the real parts of V(s) of these two models also 
agree very well with its Taylor's series expansion, using its first 4 terms (not shown). 

The results of the imaginary parts of the one and two-loop ChPT calculations are also 
shown in Fig. 1; it is seen that they differ significantly from the data. At low energy, 
ChPT results for the real parts are not bad as can be seen from Fig. 2. This is due to the 
dominance of the real subtraction constants (uncalculable in ChPT scheme) and cannot 
be used to support the validity of the ChPT as will be discussed below. 

At one-loop level, ChPT cannot be used to calculate the rms radius because the loop 
integral for this quantity is divergent. Using the measured r.m.s. as an input, the one 
loop ChPT for the pion form factor is given by Eq. ([I]). The coefficient cf hPT , where 
the subscript refers to the one-loop level is cf hPT = (9607r 2 /^m^.) _1 ~ 0.626GeV~ A and 
d ChPT = (i34407r 2 /2 m 4)-i _ 2.30GeV~ 6 The value for c calculated by the one-loop 
ChPT is in agreement with a previous determination [|18[] and is too small by a factor of 
7 compared with the sum rule value Eq. (|7|); d is also a factor of 4 too small compared 
with the sum rule value, Eq. (0). 

At two- loop level [EL [| , ChPT cannot be used to calculate the expression for c because 
of the degree of divergence of the loop integral. One can calculate, however, the expression 
for d. Although an analytical formula can be given, it is simpler to give it numerically: 
(jChPT _ ^ \Q e y-® which is a factor of 2.5 too small compared with the sum rule value, 
Eq. (0) . The disagreement with the sum rule value is now less than that from the one- loop 
calculation. 

In the following, one follows the standard dispersion relation analysis of the data, i.e. 
given the imaginary part of an amplitude, one can calculate its real part by dispersion 
relation. For this purpose, it is important to realize that the equivalent use of the Taylor's 
series is the dispersion integral (DI) on the R.H.S. of Eq. ([|). Because of the use of the 
elastic unitarity, this analysis is at best valid to a maximum of energy of 0.8 — 0.9GeV . 

One should compare the experimental values of the real part of the DI with those given 
by models, because they are the direct results of the calculation schemes, unmasked by 
the dominance of the subtraction constants at low energy. The real parts of the DI can be 
calculated using ImV(s) from experimental data, models 1, 2 and ChPT results. In Fig. 
3, for clarity, only the two-loop ChPT calculation is plotted. It is seen that the two-loop 
ChPT results, are too small compared with the corresponding real parts calculated from 
the data, not only for small s but also for large s (using n = 3 for the DI). At very low 
energy one recovers the results of the Taylor's series as discussed above. 

The real parts of the DI, calculated using models 1 and 2, are in good agreement with 
the data at low energy, but show a small deviation from the data above 0.6GeV. This 
result is expected because models 1 and 2 violate the phase theorem by a small amount 
above 0.9GeV, as mentioned above. 

It is seen that perturbation theory is inadequate for the vector pion form factor cal- 
culation for an energy below the p mass and even at the point s = which are fairly 
far away from the p mass. It would be useful to ask how large the hypothetical p mass 
would be in order that perturbation theory could be trusted. It is difficult to answer this 
question in general but one can take the large Nf model [f2D| as a guide, where Nf is the 
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number of flavors. In this model, the explicit solution for the form factor as well as for 
the scattering amplitude can be given. The expression for the form factor is given by Eq. 
(|3|). In order that perturbation theory to be valid, it is required that the coefficient c 
calculated by the ChPT to be much larger than that given by the large Nf model. This 
condition yields: 

s R » VoWfam* ~ 1.26GW 2 (9) 

or roughly s p >> 1.26GeV 2 , a condition cannot be satisfied by the physical value of the 
p mass. 

The situation may not be as bad for the scalar form factor as can be seen by considering 
also the large Nf model. Instead of Eq. @ we now have a condition on the scalar r.m.s. 
of the pion, 1/6 < r 2 scalar >= l/s% alar : 



s scaiar >> ^1920/197^^- ~ OAlGeV 2 (10) 

Using the "experimental" value of the scalar r.m.s. radius, < r 2 >= 0.6 fm 2 [|], one 
has 4 m ' ar ~ OAOGeV 2 which is equal to the R.H.S. of Eq. (0), instead of being much 
larger. The situation is better here than the vector case. It might be barely possible to 
apply the ChPT for the scalar form factor but the accuracy of the perturbative approach 
could be questionable. 

In conclusion, because of the inedequacy of the perturbative loop calculations, one 
could try two different approaches to this problem. A more radical approach consists 
in considering effective lagrangians just as it was invented for, i.e. low energy theorems; 
instead of doing perturbative loop calculations, one could try to analytically continue these 
low energy theorems to the time-like region with the technique of dispersion relation and 
the constraint of unitarity. A more standard approach is to do perturbative calculation 
but resumming the perturbative series in order to satisfy unitarity. In both approaches, 
one could use for example the inverse amplitude, N/D and Pade methods |7|, |J. This last 



method p2| , based on the results of the perturbation calculation and briefly reviewed here, 
could be an useful method but could also give a wrong answer in a number of problems if 
care was not taken. In general, one does not have a unique prescription of how to handle 
the strong interaction problem. 

The author would like to thank Dr. Torben Hannah for a detailed explanation of 
his calculation of the two-loop vector pion form factor and also for a discussion of the 
experimental situation on the pion form factor data. Useful conversations with T. N. 
Pham are acknowleledged. 
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Figure Captions 

Fig.l : The imaginary part of the vector pion form factor ImV(s), given by Eq. @, 
as a function of the energy in the GeV unit. The experimental data of the modulus of the 
form factor and the P-wave tctt phase shifts are taken from the refences [11, 12, 13. 14. |15|| . 
The solid curve is the the experimental results with experimental errors; the long-dashed 
curve is the two-loop ChPT calculation, the medium long-dashed curve is the one-loop 
ChPT calculation, the short-dashed curve is from the modified unitarized one-loop ChPT 
calculation, Eq. (§), and the dotted curve is the unitarized two- loop calculation of Hannah 





Fig. 2 : The real parts of the pion form factor as a function of energy. The curves are 
as in Fig. 1. 

Fig. 3 : The real parts of the dispersion integral ReDI as a function of energy (using 
n = 3) in Eq. The curves are as in Fig. 1; the ChPT one- loop result is not shown. 
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